Dilute limit of a strongly-interacting model of spinless fermions and hardcore bosons 

on the square lattice 
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In our model, spinless fermions (or hardcore bosons) on a square lattice hop to nearest neighbor 
sites, and also experience a hard-core repulsion at the nearest neighbor separation. This is the 
simplest model of correlated electrons and is more tractable for exact diagonalization than the Hub- 
bard model. We study systematically the dilute limit of this model by a combination of analytical 
and several numerical approaches: the two-particle problem using lattice Green functions and the 
t-matrix, the few-fermion problem using a modified t-matrix (demonstrating that the interaction 
energy is well captured by pairwise terms), and for bosons the fitting of the energy as a function 
of density to Schick's analytical result for dilute hard disks. We present the first systematic study 
for a strongly-interacting lattice model of the t-matrix, which appears as the central object in older 
theories of the existence of a two-dimensional Fermi liquid for dilute fermions with strong interac- 
tions. For our model, we can (Lanczos) diagonalize the 7x7 system at all fillings and the 20 x 20 
system with four particles, thus going far beyond previous diagonalization works on the Hubbard 
model. 

PACS numbers: 71.10.Fd, 71.10.Pm, 05.30.Jp, 74.20.Mn 
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A. The spinless fermion model 
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Since the discovery of the high-temperature supercon- 
ductors in 1986, there has been intense study of a number 
of two-dimensional models that are believed to model the 
electronic properties of the Cu02 plane of the cuprate su- 
perconductors, for example, the Hubhard model, the J 
model, and the Heisenberg model. H'El Two-dimensional 
quantum models with short-range kinetic and interac- 
tion terms are difficult to study. In one dimension, there 
are exact solutions using the |Bethe ansatz and a host of 
related analytical techniques,™ and there is a very accu- 
rate numerical method, the density-matrix renormaliza- 
tion group (DMRG),a that can be applied to large sys- 
tems relatively easily. In two dimensions, on the other 
hand, there are few exact solutions (one famous nontriv- 
ial case is the Hubbard model with one hole in a half-filled 
background, the Nagaoka stateQ), and current numerical 
methods are not satisfactory (quantum Monte Carlo is 
plagued by the negative sign problerru at low tempera- 
tures and at majiy fillings of interest and the DMRG in 
two dimensionsEl is still in early development stage). 

The most reliable method for studying complicated 
quantum systems is exact diagonalization, which means 
enumerating all basis states and diagonalizing the result- 
ing Hamiltonian matrix. Of course, this method is com- 
putationally limited by the growth of the Hilbert space 
which is in general exponential in the number of parti- 
cles and the lattice size. The 4x4 Hubbard model with 
16 electrons, 8 spin-up and 8 spin-down, after reduction 
by particle conservation, translation, and the symmetries 
of the_square, has 1,310,242 states in the largest matrix 
block,Ll and can be diagonalized using the well-known 
Lanczos method. The Hubbard model has been diag- 
onalized for the 4x4 lattice (see. e.g., Ref. ||), and at low 
filling (four electrons) for 6 x 6lI with extensive employ- 
ment of symmetries. 



We have asked the question: Is there a model that 
contains the basic ingredient of short-range hopping and 
interaction hut is simpler, in the exact diagonalization 
sense, than the Hubbard model? The answer is yes: we 
can neglect the spin. We obtain the following Hamilto- 
nian for spinless fermions. 



H^-t 



4c j 



(1.1) 



where c\ and Ci are spinless fermion creation and anni- 
hilation operators at site i, fii — cjci the number oper- 
ator, t the nearest-neighbor hopping amplitude, and V 
the nearest-neighbor interaction. Note that with spinless 
fermions, there can be at the most one particle per site; 
no on-site interaction (as that in the Hubbard model) 
is possible, and we have included in our Hamiltonian 
nearest-neighbor repulsion. 



The spinless fermion model, Eq. (1.1), is a two-state 
model, and the number of basis states for a iV-site system 
is 2^, which is a significant reduction from the 4''^ of the 
Hubbard model. We can further reduce the number of 
basis states by taking the nearest-neighbor interaction 
V = -\-oo, i.e., infinite repulsion, which excludes nearest 
neighbors, giving roughly 2^/^ states. 

Th e s pinless fermion model with infinite repulsion 
Eq. (^]^) contains a significant reduction of the Hilbert 
space. After using particle conservation and translation 
symmetry (but not point group symmetry), the largest 
matrix for the 7x7 system has 1,906,532 states (for 
11 particles), and we can therefore compute for all fill- 
ings the 7x7 system whereas for the Hubbard model 
4 X 4 is basically the limit. This of course means that for 
certain limits we can also go much further than the Hub- 
bard model, for example, we can handle four particles 
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on a 20 X 20 lattice where the number of basis states is 
2, 472, 147. This extended capability with our model has 
enabled us to obtain a number of results that are difhcult 
to obtain with the Hubbard model. 

An added feature of our model is that the basis set 
for the spinless fermion problem is identical to that for 
the hardcore boson problem, because with hardcore re- 
pulsion, there can be at the most one boson at one site 
also. Therefore, without computational difRculty, we can 
study numerically both the spinless fermion and hardcore 
boson problem. 

Spinless fermions can also be realized in experiments, 
for example, the spin polarized '^He due to a strong mag- 
netic field, or ferro or ferri-magnetic electronic systems 
where one spin-band is filled. The one-dimensional spin- 
less fermion model with finite repulsion is solved exactly 
using Bethe ansatz.t2l The infinite-dimensional problem 
is studied in Ref. ^ A very different approach using the 
renormalization group for fermions is done in Ref. 12, A 



Monte Carlo study of the two-dimensional model at half- 
filling only and low temperatures is in Ref. ^ which, 
dating back to 1985, is one of the earliest quantum Monte 
Carlo simulations for fermions. (It is no coincidence that 
they chose the model with the smallest Hilbert space.) 

Considering the tremendous effort that has been de- 
voted to the Hubbard model and the close resemblance 



of our model, Eq. (1.1), to the Hubbard model, it is sur- 
prising that works on this spinless model have been rather 
sparse, though it has been commented tiiat the spinless 
model offers considerable simplifications liil 

This paper is one of the two that we are publish- 
ing to study systematically the two-dimensional spinless 
fermion and hardcore boson model with infinite nearest- 
neighbor repulsion. The present paper focuses on the 
dilute limit, troating the problem of a few particles, and 
the otiier paperllJ will focus on the dense limit, near half- 
filled,Ej where stripes (that are holes lining up across the 
lattice) are natural objects (see Ref. |l^ for a condensed 
study of stripes in this model). We will use Lanczos ex- 
act diagonalization, exploiting the much-reduced Hilbert 
space of our model, and a number of analytical tech- 
niques, for example, in this paper, lattice Green func- 
tions and the t-matrix. One of the goals of these two 
papers is to advertise this model of spinless fermions to 
the strongly-correlated electron community, as we believe 
that it is the simplest model of correlated fermions and 
deserves more research effort and better understanding. 

The prior work most comparable to ours may be the 
studies of four spinless electrons in a_£ x 6 lattice, with 
Coulomb repulsion, by Pichard et aljtB their motivation 
was the Wigner crystal melting and the competition of 
Coulomb interactions with Anderson localization when a 
disorder potential is turned on. 



B. The t-matrix 

At the dilute limit of our model, the scattering t-matrix 
is of fundamental importance. For two particles, we ex- 
pect that, at least when the potential V is small, we can 
write a perturbative equation for energy. 



(1.2) 



which is to say that the exact interacting energy of two 
particles is the noninteracting energy £{q_i) +£{q2), for 
a pair of momenta qi and q2, plus a correction term 
AE due to the interaction V. And with more than two 
particles, at least when the particle density is low, we 
expect to have 



£; = 5^£(q) + i^Ai?(q,q'). 



(1.3) 



Eq. ( |l.3| ) is central in Fermi liquid theory, where it is 
justified by the so-called "adiabatic continuation" idea, 
which says that interacting fermion states correspond 
one-to-one to noninteracting ones as we slowly switch on 
a potential. 

In the boson case, because many bosons can occupy 
one quantum mechanical state and form a condensate, 
Eq. (1.3) shou ld b e modified, but with only two bosons, 
we expect Eq. (1^) should be valid (in that t he c orrection 
vanishes in the dilute limit). Eqs. (1.2) and (1.3) are used 
when we look at a list of noninteracting energies and draw 
correspondences with the interacting energies, the energy 
shift being packaged in the term AE. 

One p ossible ob jection to the above formulas 
(Eqs. ( |l.2| ) and ( |l.3| )) is that they appear to be pertur- 
bative, yet the interaction potential in our problem is 
infinitely strong, so the first-order (first Born approxi- 
mation) scattering amplitude, being proportional to the 
potential, is infinite too. However, this singular poten- 
tial scattering problem (e.g., hard-sphere interaction in 
3D) has been solved (see Ref. by replacing the po- 
tential with the so-called scattering length, which is fi- 
nite even when the potential is infinite. As we review 
in Appendix a perturbation series (Born series) can 
be written down (that corresponds to a series of the so- 
called ladder diagrams) and even though each term is 
proportional to the potential, t he s um of all terms (the 
t-matrix, AE in Eqs. ([ij) and (U)) is finite. 

Because the t-matrix captures two-body interaction ef- 
fects, it is the centerpiece of dilute fermion and boson 
calculations with strong interactions. Field-theoretical 
calculations in both three and two dimensions are based 
on the ladder diagrams and the t-matcipc. See Fetter and 
WaleckaEJ for the 3Iljproblem, Schickcil for the 2D boson 
problem and Bloomcj the 2D fermion problem. For lat- 
tice fermion problems, Kanamoric^ derived the t-matrix 
for a tight-binding model that is essentially|-a Hubbard 
model (this work is also described in YosidacJ). And in 
Ref. |2^, the t-matrix is worked out explicitly for the Hub- 
bard model, and Kanamori's result is obtained. Ref. 
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also evaluated the t-matrix for the dilute limit in three 
dimensions and obtained a functional dependence on par- 
ticle density. j— , 

Rudin and MattiaSj used the t-matrix expression de- 
rived in Refs. and ^ and found upper and lower 
bounds of the fermion t-matrix in two dimensions in 
terms of particle density. Rudin and Mattis's result for 
the low-density limit of the two-dimensional Hubbard 
model is of the same functional form as Bloom's dia- 
grammatical, calculation for the two-dimensional fermion 
hard disks.E3 Since the discovery of high-temperature su- 
perconductors, Bloom's calculation has received a lot of 
attention because of its relevance to the validity of the 
Fermi liquid description of dilute fermions in two dimen- 
sions. There have .bejen a number of works on the 2D di- 
lute Fermi gaspZlSo and on the dilute limit of 2D Hub- 
bard modelJ23 all using the t-matrix, but these results 
have not been checked by numerical calculations. 

In fact, we are not aware of a systematic study of the 
t-matrix for a lattice model. In this paper, we present the 
first such study for the two-particle problem in Sec. [II 
(for bosons and fermions) and the few-fermion problem 
in Sec. IV. We check the t-matrix results with exact 



diagonalization data and show that our t-matrix on a 
lattice is the sum of the two-body scattering terms to 
infinite order. 



C. Paper organization 

In this paper, we will study systematically the dilute 
limit of our model Eq. (O), focusing on the problem of 
a few particles. Our paper is divided into four parts. 

In Sec. H, the two-particle (boson and fermion) 
problem is studied. We formulated the two-particle 
Schrodinger equation using lattice Green functions, em- 
ploy some of its recursion relations to simplify the prob- 
lem, and obtain the two-boson ground state energy in 
the large-lattice limit. Using the two-particle result, we 
then study the problem of a few particles and obtain an 
expression for ground state energy on a large lattice. 

In Sec. Ill, the two-particle problem is then cast into a 
different form, emphasizing the scatterings between the 
two particles. The result is the t-matrix, that is ex- 
act for the two-particle problem and contains all two- 
body scattering terms. We will study the two-particle 
t-matrix in great detail, showing the differences between 
the boson and fermion cases, and demonstrating that the 
first t-matrix iteration is often a good approximation for 
fermion energy. In Appendix we show explicitly that 
the t-matrix we obtain is the sum total of all two-body 
scattering terms. 



The problem of a few fermions is taken up in Sec. IV 



First, the fermion shell effect is discussed and demon- 
strated from diagonalization, and we show the difference 
for bosons and fermions. We show the modifications to 
the two-fermion t-matrix that enable us to calculate en- 
ergies for three, four, and five particles. Using this t- 



matrix, we can compute the interaction corrections to 
the noninteracting energy and trace the change in the 
energy spectrum from the nointeracting one to the inter- 
acting one. 

Finally, in Sec. 0, we discuss the energy per parti- 
cle curve for dilute bosons and fermions. We have|-stud- 
ied the two-dimensional results derived by SchiclH for 
bosons and Bloomcj for fermions by fitting the data from 
diagonalization for a number of lattices. Schick's result 
for dilute bosons is checked nicely, and we explain that 
for spinless fermions in our model we will need the p- 
wave scattering term, which is not included in Bloom's 
calculation. 

In Appendix |^ we discuss briefly our exact diagonal- 
ization computer program, which can handle arbitary pe- 
riodic boundaries specified by two vectors on the square 
lattice and uses translation symmetry to reduce the ma- 
trix size. 



II. THE TWO-PARTICLE PROBLEM 

The two-particle problem has appeared in many dif- 
ferent contexts. The most familiar one is the hydro- 
gen atom problem in introductory quantum mechanics 
textbooks. The two-magnon problem is closely related 
mathematically to our two-particle problem, and it has 
been solved in arbitary dimensions for ferromagnets (see 
e.g., Ref. Another important two-particle prob- 

lem is the Cooper problem, with two electrons in the 
presence of a Fermi sea (see e.g., Ref. |3^). And mo- 
tivated by the possibility of Cooper pair formation in 
high-temperature superconductors, there have also been 
a number, nf stjjjliejji on bound states on a two-dimensional 
lattice.HSMy'ElEl The two-electron problem in the 
plain two-dimensional repulsive Hubbard model is stud- 
ied in Ref. and ground state energy in the large-lattice 
limit is obtained analytically. 

In this section, we present a rather complete calcula- 
tion for the two-particle problem in our model, treating 
both bosons and fermions. With infinite repulsive in- 
teraction in our model, we are not interested in finding 
bound states. We calculate eigenenergies for all states for 
a finite-size lattice, and our caiculation is more compli- 
cated than the Hubbard modeE3 case because of nearest- 
neighbor (in place of on-site) interaction. Where the 
Green function in the Hubbard case was a scalar object, 
in our case it is replaced by a 4 x 4 matrix, correspond- 
ing to the four nearest neighbor sites where the potential 
acts. This Green function study of the two-particle prob- 
lem is closely related to the treatment of the two-electron 
problem in the Hiibbard modelEj and that in an extended 
Hubbard model.E3 We will show the use of lattice symme- 
try and recursion relations to simplify the problem with 
nearest- neighbor interactions. 



4 



A. Preliminary 



In this two-particle calculation, wc will work in mo- 

-start with a Hamiltonian 



mentum space, and we wil 
more general than Eq. (p, . 1|) , 

H = T + U, (2.1) 
T = ^t(r2-ri)cl;^cr,, (2.2) 

C/ = E^(^2-ri)ct^Cr,ct^Cr,. (2.3) 

rir2 

Here we have allowed hopping and interaction between 
any two lattice sites, but we require that both depend 
only on the separation between the two vectors and both 
have inversion symmetry. That is t(ri,r2) ~ t{r2 — ri), 
t(-r) - t{r), F(ri,r2) = Vill-ri), and Vi^r) = y(r). 
In momentum space, Eqs. (2.2) and (|2.3|) become. 



T = 



U 



^£(p)4cp, 
p 

\j J2 ^(k)c]34,Cp/+kCp_k, 



2N 



pp'k 



where 



f(p)=^^(r)e^P^ 

r 



(2.4) 
(2.5) 

(2.6) 
(2.7) 



with ^(-p) = £(p) and V(-k) = V(k). Eqs. (|2j) 
and (|2.3| ) reduce to our nearest-neighbor Hamiltonian 
Eq. (|1 . 1|) if we take, 



t{r)^\-'' r = (±l,0)(0,±l), (2.8) 
0, otherwise. 



V{r) 



_jV, r-(±l,0)(0,±l), 
0, otherwise. 



(2.9) 



where we have taken the lattice constant to be unity, 
and the nearest-neighbor vectors will be called 



Ri = (1,0), R2 = (-1,0), R3 
Then we have, 



(0,1),R4 = (0,-1). 

(2.10) 



£(p) = —2t{cospx + cospy), (2-11) 
V{k) = 2V{coskx+ cos ky). (2.12) 

Note that the structure of later equations depends sensi- 
tively on having four sites in Eq. (2^) where V{ r) 7 ^ Om, 
but does not depend much o n the form of Eq. (2.8) and 
the resulting dispersion Eq. (2.11). 



Using momentum conservation of Eq. (2.1), the two- 
particle wave function that wc will use is. 



= ^g(q)|q,P-q), 



(2.13) 



where the sum is over the whole Brillouin zone, and the 
coefficient ^(q) satisfies, 

.g(P - q) - Sbfgiq), (2.14) 

where Sbf = ^ for bosons and —1 for fcrmions. 

B. Green function equations 

Applying t he m ore general form of the Hamiltonian 
operator Eqs ( ^ ) and ( ^^ to the state Eq. ( |2.13| ), the 
Schrodinger equation {E — T)\ip) = U\ip) becomes 

(£;-£(q)-£:(P-q))5(q) = ^ J] ^(q-k).9(k). (2.15) 



Eq. ( |2.15| ) is a matrix equation Ag = Eq where Aqk = 
(£(q) +£{P- q))^qk + V{q - k)/N. If V is not infinity, 
this N X N matrix A can be diagonalized, and E and 
f;(q) are respectively the eigenvalue and eigenvector. To 
deal with V = +00, we need some further manipulations. 

We consider the case when E ^ f(q) + i^^(P — q), for 
any q, which is to say, the energy E is not the energy of 
a noninteracting pair. The (lattice) Fourier transform of 
the coefficients g{q) is 



g(r) = ^e-''-g(q); 



(2.16) 



this is just the real-space wavefunction in terms of the 
relative coordinate r. Define the lattice Green function. 



1 --^ p«q-(r'-r) 

G(^,P.-.--) = ^E ^_,(,)_,(P_,) 



(2.17) 



then after dividing the first fact or from both sides of our 
Schrodinger equation, Eq. ( 2.15 ), and Fourier transform- 
ing, we obtain 
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(r) = ^G(i?,P;r,r')nr').9(r')- 



(2.18) 



In the following we ret urn to the nearest-neighbor po- 
tenti al V' (r) in Eq. ( ^.9[ ). The Green function sum in 
Eq. (|2.18| ) then has only four terms. 



~g{r) = J2G{E,P;r,R,)iVg{R,)), 



(2.19) 



summed over the separations in Eq. ( 2.10 ). If 1*^ also 
restr ict r to the four nearest-neighbor vectors J13 then 
Eq. ( |2.19D becomes, 
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(2.20) 



If we define the 4x4 matrix, 

^„(S,P) = G(S,P;R„Rj), (2.21) 
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and a 4 X 1 vector 
matrix equation, 



g(Rj), then we obtain a simple 



{I -g{E,P)V)(f)^0. 



(2.22) 



We can also rewrite this equation as an equation for en- 
ergy using the determinant, 

det{I -g{E,P)V) ^0. (2.23) 

With V = +CXD, we have even simpler equations 

g{E,P){V^) ^0, (2.24) 

and 

dete;(£;,P) = 0. (2.25) 

Notice we write V(p to denote the lim it as V oo; i t 
would not do to write simply cj) in Eqs. ( 2.22 ) and ( ^.24 ), 
since ^ as ^ oo (being the amplitude of the rela- 
tive wavefunction at the forbidden separations {R^}.). 

For the Hubbard model, there is only on-site interac- 
tion, so V{r ) is nonzero only when r = 0, and the sum 
in Eq. (2.18) has only one term. Eq. ( 2.2C| ) is simply a 
scalar equation, which, after q can cels from both sides of 
the equation and using Eq. ( 2.17 ), gives, 



V 

N 



E 



E - Siq) - £{P - q)' 



(2.26) 



which is exactly the result in Ref. p7 



1. Simplifications for rectangular boundaries 



We specialize to the case of total momentum P = and 



recta ngula r-boundary lattices. We have from Eqs. ( |2.17] ) 
and d^TI), 



G^.iE) 



1 ^ COs(ga(^jx 



Rix)) COa{qy{Rjy - Riy)) 



N 



E - 2£{ci) 

(2.27) 

where t he p otential is nonzero on the sites {Ri} given 
by Eq. ( 2.1C ) and in the last step we have used the sym- 
metry properties of the dispersion relation £ (q^ , Qy ) = 
£iqx,-qy) = £{-qx-,qy)- Obviously Eq. (|2.27| ) is a 



function of displacements Rj — Ri, which (in view of 
Eq. ( |2.1C| ) can be (0,0), (1,1), (2,0), or any vector related 
by square symmetry. It is convenient for this and later 
sections to define a new notation for the Green function 
Qij, emphasizing its dependence on Rj — R^ = (m,n), 



r{E,m,n) ^ 



cos(mq^) cos{nqy) 
N ^ E + 4: cos qx + 4: cos qy ' 



(2.28) 



where the sum is over the N wavevectors q = 
{2T:Lj;/Lx,2'Kly/Ly) with < Ix < and < ly < Ly 
(for one Brillouin zone), and we have used the expression 
for £(q) from Eq. (2.11) (and taken t = 1). 



This Green function for rectangular-boundary lattices 
satisfies the following reflection properties, 

r(i5, m, n) — T(E, —m, n) 

= T{E,m,~n) ^T{E,-m,-n). (2.29) 

And if we have a square lattice (L^ ^ Ly) we also have 

r{E,m,n) ^T{E,n,m), (2.30) 

Eq. ( 2.27| ) can be written as, 



Giji^j) — r(-E, Rjx — Rix,Rj 



Ri' 



(2.31) 



Us ing the reflection properties of T{E,m,n)^ 
Eq. ( |2.29| ), and the definition Eq. ( |2.3lD , our Green 
function matrix becomes, 



/ a c b b 

c a b b 

b b a d 

\ b b d a 



(2.32) 



where a = r{E, 0, 0), 6 = T{E, 1, 1), c = T{E, 2, 0), and 
d — T{E, 0,2). The eigenvalues and eigenvectors of this 
matrix are, 

Xfi=a-c, (1,-1,0,0) 
Xf2=a-d, (0,0,1,-1) 



c + d ^1662 + (c-d)2 
A6i,62 — a H ^ — ± ^ : 

(Wl,2,l'l,2, 1, 1) 



(2.33) 



where vi and V2 are complicated functions of a, b, c, and 
d. 

The exact energy E makes the matrix Qij [E) sin- 
gular, which means that one of the eigenvalues has to 
be zero. From Eq. ( ^.24 ), the null eigenvector of G is 



= F(g(^,5(R2),g(R3),5(R4)), in terms of {RJ 
as in Eq. ( ^.1C| ). The relative wavefunction should be 
odd or even under inversion , depending on statistics, 
i.e. g{ —r) = Sfcf e' ^ '".g(r ) which follows immediately from 
Eqs. (^.16 ) and ( 2.14 ). Inversion, acting on nearest- 
neighbor vectors Eq. ( 2.10| ), induces Ri ^ R2 and 
R3 ^ R4; thus with P = 0, we should have V(j)i = —V^ii 
and V^03 — —V(^^ for fermions, and V(\)\ — V(j)2 and 
V(p3 — V(t)i for bosons. Inspecting the eigenvectors we 
obtained in Eq. (2.33), we see that those corresponding to 
A/1 j2 are antisymmetric under inversion - corresponding 
to a "p- wave- like" (relative angular momentum 1) state 
for fermions. Setting A/i =0, we get a = c, or setting 
A/2 = a = d, which respectively mean 

r(£;,o,o) -r(i;,2,o) = o,or (2.34) 
r(£;,o,o) -r(i;,o,2) = 0. (2.35) 

Associated with the even eigenvectors are Abi_;,2 which 
are identified as boson eigenvalues. 
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2. Simplifications for square boundaries 



The boson eigenvalues, Eq. (2.33), are rather comph- 
cated for general rectangular-boundary lattices. For a 
square-boundary lattice, using Eq. (2.30), we get c = d 
in the matrix Eq. ( ^.32 ), which makes the fermion eigen- 
value s Xf iF2 degenerate. The boson eigenvalues in 
Eq. (2.33) simplify greatly to Afci = a + 2b + c and 
Xb2 = a — 26 + c, which means that the boson energy 
equations are, 



r(£;,o,o) 

T{E, 0,0) 



2r(i?,i,i) 



r{E,2,0) = 0, (2.36) 
r(£;,2,0) = 0. (2.37) 



The corresponding eigenvectors simphfy too, to (1, 1,1,1) 
and (1,1,— 1,-1) respectively, which may be described 
as "s-wave-like" and "d-wave-like" , i.e. relative angular 
momentum and 2. 



C. Large-L asymptotics for two-boson energy 



Eqs. ( |2.34| ), ( |2.35| ), ( |2.36| ) and ( |2.37| ) are much better 
starting points for analytical c alcul ations than the origi- 
nal determinant equation Eq. (2.25). In the center of the 
problem is the lattice Green function T{E,m,n) defined 
in Eq. (2.28). Many of the lattice calculatiopspjmp-dfljp. 
to evaluating these lattice Green functions .OoE3EjE3'EJ 
In this section, we derive the large-lattice two-boson en- 
ergy using the recursion and symmetry relations of the 
Green function T{E, m, n). 

The Green function T{E, m, n) for general m and n and 
finite lattice are difficult to evaluate. The good thing is 
that there are a number of recursion relatinijis ponnecting 
the Green functions at different m and n£3o These are 
trivial to derive after noting that Eq. (2.17) (for P = 0) 
can be written 

[^-(4+A2)-(4+A;2)]G(S,0;r,r') = d,=oSr'=o (2.38) 

where is the discrete Laplacian, (A^ -I- 4)/(r) = 
^j/(r -I- Ri) for any function /(r), where the sum is 
over neighbor vectors Eq. (2.10). The two recursion re- 
lations that we will use are 

£;r(£:,o,o) + 4r(£;, i,o) + 4r(£;,o, 1) = 1,(2.39) 

r{E, 0,0) + 2 T{E, 1, 1) + T{E, 2, 0) 

+ ^ETiE,l,0) = 0.(2.40) 



Usin g Eqs . ( |2.30| ), ( p.39D , and ( |2.40D , the boson equation 
Eq. (2.36) for square-boundary lattices simplifies to 



T{E, 0,0) 



E' 



(2.41) 



with eigenvector (1,1,1,1). 

Next we compute the leading form of T{E, 0,0) for 
large L of a square-boundary lattice. The calculation is 



close to that in Ref. ^ for the Hubbard model. We define 
E = —8 -I- AE. Because the lowest energy of an indepen- 
dent particle is £{0) = —4, AE' is the energy correction 
to two independent parti cle e nergy at zero momentum. 
Then we have, from Eq. (2.28), 



r(£;,o,o) 



-E- 



1 



iV ^ + 4 COS -f 4 cos qy 



47V ^ 2 - cos qx - COS qy - AE/A ' 
1 



]_ r dq 

1 InL 

const. 



L'^AE Att 



(2.42) 



We should discuss the number of approximations we have 
made to extract this leading dependence in L. First ex- 
cept in the q = term we have ignored the AE term, 
assuming it is small as compared to (with q 7^ 0). This 
is justified as we only want the leading term in the large- 
L limit. Using an integral for a lattice sum is another 
approximation. We choose the lower limit of integration 
to be 2tt / L corresponding to the first wavevectors after 
(0, 0) is taken out of the sum. We also used the quadratic 
approximation for the energy dispersion f (q) appearing 
in the denominator. 

Using the boson energy equation Eq. (^.41 ) and the 
large- 1/ limit of the Green function Eq. (2.42), we get. 



1 



-8 + A£; L^AE 



InL 
47r 



const. 



(2.43) 



In the large-L limit, AE (as it is the interaction 
correction to the noninteracting energy), so we get, to 
the leading order of L, 



AE = 



4:77 



L^lnL 



(2.44) 



We will check Eq. {^A^ in Sec. pD 



D. Large-L asymptotics for few- particle energy 



The procedure used in Sec. p C] for two bosons can also 
be applied to problems with a few particles. For a few 
particles on a large lattice with short-range (here nearest- 
neighbor) interaction, two-particle interaction is the main 
contribution to energy. We write for two particles. 



E{2,L) = Eo{2,L) + AE{L). 



(2.45) 



Here in this section we use the notation E{M, L) and 
Eq{M,L) to denote the M-particle exact and noninter- 
acting ground state energies respectively and emphasize 
the dependence of AE on L by using AE{L). It is rea- 
sonable to expect that the energy for M particles is the 
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noninteracting energy plus interaction corrections from 
the M(M — l)/2 pairs of particles. We then have, 



E{M, L)^E„{M, L) 



M{M - 1) 



AE{L). (2.46) 



For bosons, Eq{M,L) = —AM, because in the ground 
state, all bosons occupy the zero-momentum state. On 
the other hand, for fermions, because of Pauli exclusion, 
no two fermions can occupy the same state, the non- 
interacting ground state is obtained from filling the M 
fermions from the lowest state (k — 0) up. 

Eq. (^.46| ) imphes that plotting 2{E{M,L) - 
Eo{M, L))/{M{M - 1)) versus L for different M should 
all asymptotically at large L approach Ai?(L). In Fig. ||, 
we do such plots, for bosons and fermions with M = 
2, 3, 4, 5. The fermion results, from p-wave scattering (as 
our spinless fermion wave function has to be antisym- 
metric), are much smaller than the boson results (bold 
curves) from s-wave scattering. 
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FIG. 1; Boson andfermion2(£(Af,L)-So(M,L))/(M(Af- 
1)) versus L for AI — 2, 3, 4, 5. All curves appear to converge 
at large L. The fermion (p-wave) result is much less than the 
boson result (s-wave). The M = 4 plot goes to L = 20 and 
the M = 5 plot to L = 10. The boson M = 5 curve is too 
high to be included in this plot. 



Note that in our calculation for r{E,0,0) Eq. ( ^.42[ ), 
we have neglected the contribution of AE in the denom- 
inator except for the first term (q = 0). Now with the 
leading form of AE E n. (^.441) , we can obviously plug 
E -8 + AE into Eq. j2.42| ) to get the form of the next 
term. 



AE 



An 



InL 



A 



C 



B 



InL (InL) 2 



(2.47) 



A + B- 



InL 



C 



1 

h^ 



(L;o(M,L) = -4M), 

{E{M, L) + AM)L^ In L 
2-kM{M - 1) 

(2.48) 

In Fig. |, we plot {E{M, L)+AM)L^ \n L / {2tt M {M - I)) 
versus 1/ InL for M = 2, 3, 4, 5, using the boson data in 
Fig. [l]. Quadratic polynomial fitting is done for M = 2, 3, 
where we have more data than M ~ 4, 5. The c oeffi- 
cient A « 1 for both fits, implying, from Eq. (2.47), the 
leading-order term in AE{L) is indeed Att/ (L^ InL). B 
and C from two fits are also comparable. 



3.5 



3.0 



2.5 



2.0 



1.5 



1.0 



+ M = 2, 1 .009+1 .029 X+2.628 X*X 
) M = 3, 1 .029+0.771 X+3.858 X*X 
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FIG. 2: Boson (.B(M, L) -|-4Af)L2 In L/(27rM(M - 1)) versus 
1/ In L for M = 2,3,4,5. Quadratic polynomial fitting is done 
for M = 2 and Af = 3. The fitted constant coeflicients are 
approximately one, and the other coeflicients from M = 2 
and Af = 3 are comparable. 



To summarize, from Eqs. ( 2.46 ) and ( 2.47 ) and fitting 
in Fig. 1^, we find that in our model the energy of a small 
number M of bosons on a large L x L lattice is to the 
leading order of L, 



E{M, L) « -AM - 



M{M - 1) 47r 
2 L2 In L ' 



(2.49) 



For two fermions on a large L x L lattice, the noninter- 
acting energy ~ the lead term in Eq. (2.49) - is obviously 
lower for P = (0,1) than for P — (0, 0). We have not 
worked out the asymptotic behavior for P 7^ (0, 0). 



III. THE TWO-PARTICLE T-MATRIX 



In Sec. 



Using Eqs. (2.46) and (2.47), we get, for a few bosons 



11 C| and Sec. |ID|, we studied the ground state 
energy of a few particles on a large lattice, and we showed 
that the energy of M particles can be approximated by 
summing the energy of the M{M — l)/2 pairs. In this 
section, we reformulate the equations for two particles 
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and derive a scattering matrix, the t-matrix. T he t - 
matrix gives us equations of the form Eqs. (L2) and ( |l.3| ) 
which are m ore precis e sta tements of the ideas presented 
in Sec 



II C and Sec. P D|. They apply to smaU lattices 



and to excited states. 



at the same time: the symmetric solution ^(q) = g(P—q) 
is a boson solution and the antisymmetric solution g(q) = 
—g(P — q) is a fermion solution (see Eq. ( |2.14| )). 

B. Derivation of the t-matrix 



A. Setup and symmetry 

To have an equation in the form of Eq. ( |l.2| ) , we start 
with any pair of momentum vectors qi and q2 and write 
noninteracting energy of the pair Eo = 5(qi)+5(q2) and 
total mom entu m P = qi +q2. Because our Hamiltonian, 
Eqs. (2^) and ( |2.5| ), conserves total momentum, we can 
restrict our basis states to |q, P — q). It is tempting to 
take |qi,P — qi) and |q2,P — q2) as our nonperturbed 
states, but there can be other two-particle states with 
the same total momentum P and energy Eq. 



In fact, using our energy dispersion function Eq. (2.11), 
if we write qi = (qix^qiy) and q2 = {q2x,q2y), and de- 
fine qa = {qix,q2y) and q4 = iq2x,qiy) then we have, 
qi + q2 = qs + q4 and £{qi) + £(q2) = f (qs) + f (q4)- 
We call this fact, that component exchanges in the pair 
qi and q2 gives a pair qa and q4 that have the same to- 
tal momentum and energy, the pair component exchange 
symmetry of our energy dispersion function f (q). This 
symmetry is is due to the fact that our £(q) is separable 
into a X part and a y part (i.e., f (q) = £x{qx)r-+ £y{qy) 
where Sxiq) = —2tcosq — £y{q) in our model) .E3 

The pair component exchange symmetry says that if 
qix ^ q2x and qiy ^ q2y, then the state |q3,q4), with qa 
and q4 defined above using component exchange, has the 
same total momentum and energy as |qi,q2)- The de- 
generate perturbation theory requires |q3,q4) should be 
included in the set of nonperturbed states with |qi, q2)- 

With a noninteracting two-particle energy Eq and to- 
tal momentum P, we divide the N wavevectors into two 
disjoint sets, 

Qo - {q I f(q) + f (P - q) = i^o}, - {q I q ^ Oo}. 

(3.1) 

Note that if q g Qo then P — q e Qq. Denote A^o the 
number of elements in Qq and N — N — Nq the number of 
elem ents in Q. With this separation of q, our eigenstate 
Eq. ( |2.13[) becomes 

W= 5] 5(q)|q,P-q) + Eff(q)|q,P-q), (3.2) 

qSQo qGQ 

where the first sum contains all states whose noninter- 
acting energy is degenerate. Using the idea of degener- 
ate perturbation theory, we expect to be able to find a 
secular matrix T, A'o x A'o, for the degenerate states in 
Qo only, and T will eventually be our momentum space 
t-matrix, which we will deri ve now. 

Note that using Eq. ( ^.14| ), t he n umber of independent 
states in the first sum of Eq. (3.2) is less than Nq. We 
include both |q, P — q) and |P — q, q) in our calculation 
because we are considering boson and fermion problems 



Our purpose is to derive a set of closed equation s fo r 
f;(q), the coefficent in our two-particle state Eq. ([3.2|), 
with q e Qo- 

The Schrodinger equation for the two-particle state 
lip), Eq. ( 2.15| ), can now be written as. 



(ii;-£(q)-£(P-q))g(q) =l^e''i'-V(r')g(r'), (3.3) 

r' 

wher e g(r ) is the Fourier transform of g(q) as defined in 
Eq. (|t|). _ 

For q G Q, if we assume that E ^ £{q) + £{P ~ q), 
Eq. ( |3.3|) becomes 

~g\v) = 5(r) - ^ G{E, P; r, r')l/(r')ff(r'), (3.4) 



where we have defined a Green function for the set Q, 



G{E,V-vy) 



E 



piq-(r'-r) 



TV ^i?-£(q)-£(P-q) 



and a Fourier transform with vectors in Qqj 
~g\v) = J2 e"'''"".9(q). 

qeQo 



, (3.5) 



(3.6) 



By res tric ting to the nearest-neighbor repulsion poten- 
tial Eq. (E|), Eq. §1) becomes, 



~g"{r) = ff(r) - ^ G{E, P; r, R,)V~giR,) (3.7) 



summed over neighbor vector Eq. (2.10). Now restricting 
r = Ri in Eq. (3/7), we get a set of four equations, 



^0,=<|>^-J2g^JiE,P){Vq^,), 



(3.8) 



where we have written 

(E, P) = G{E, P;-R,,-R,) (3.9) 



and = g(Ri) and — ^"(Ri). Eq. ( |3.8|) is a matrix 
equation, 



4>Q = {i-giE,p)v)^, 



(3.10) 



where ^ is 4 x 4, and 0o are 4x1, and F is a scalar 
(strength of potential) . And we can invert the matrix to 
get, 



{i-g{E,p)vy 



(3.11) 
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This is a key result, as we have expressed the desired 
function g, a Fourier transform of g(q) including all q, in 
terms of g° which includes only q S Qo; the information 
about other q € Q was packaged into the Green function 



Now we go back to Eq. (3.3), restri ct th e summation 
to Ri, and substitute in Vgi from Eq. (3.11), and we get, 



iE-£{q)~£{P~q))g{q)^ J2 T(i?, P; q, q')g(q'). 

(3.12) 

where in the la st st ep we have used the Fourier transform 
of .g"(Rj) Eq. and defined, 

r(£;,P;q,q') = 1 ^ e^qR^g-^'i'^^- {V{I ~ g{E)Vr') 
ii 

(3.13) 



If we restrict q G Qo in Eq- (3.12), then we have, 



iE-Eo)giq)^ J2 T{E,P;q,q')giq'), (3.14) 
q'eQo 

which means, 

E^ Eq + Eigenvalue(r(i;)), (3.15) 
where we have written 



%1,^'iE) ^T{E, P;q,q') 



(3.16) 



and left out the dependence on P. Tq.q' is the t- matrix 



in momentum space. Both q and q' in Eq. (3.14) are 
in Qo, which means that if there are A^o elements in Qq 
then t he ma trix T{E) is A^o x Nq. 

Eq. ( 3.15| ) is our desired equation that shows explicitly 
the interaction correction to the noninteracting energy 
Eq. In Appendix ^ we show the physical meaning of 
T{E, P; q, q') in the language of diagrammatic perturba- 
tion theory, namely it is the sum total of all the terms 
with repeated scattering of the same two particles. This 
t-matrix formalism for the two-particle problem is there- 
fore exact, and it is exactly equivalent to the Schrodinger 
equation and the Green function formalism in Sec. 
The resulting equation is an impl icit equation on E, of 
the form E = Eq + AE{E) of Eq. (|L2|), and we wiU show 
in a later section that for fermions the approximation 
E ^ Eq + AE{Eo) is often very good. 

Note also that for our case V = +oo, the t-matrix 



expression Eq. (3.13) becomes 



T(i?, P;q,q') = ^ 5^ e^'^^-e-^'i'^^ 



'GiEY 



(3.17) 

where the potential V cancels out, giving a finite value. 
This is one of the advantages of the t-matrix formal- 
ism that it can deal with infinite (singular) potential, 
for which straightforward perturbation theory would di- 
verge. 



The definition of T(£;,P;q,q') in Eq. ( |3.13| ) is a 
Fourier transform of the real space quantity V{I — 
Q{E)V)^^ . Here G is 4 x 4 because we have nearest- 
neighbor interaction. When there is only on-site interac- 
tion, as is in the usual Hu bbar d model case, G{E) — 
G(£;,P,(0,0),(0,0)), Eq. (U), is a scalar. Then, 
we can simply use the scalar quantity V/{I — GV2- 
which is-|the t-matrix that has appeared ip_Kanamori,Ej 
MattisEj Rudin and Mattis,c3 and YosidaBJ Our expres- 



sion, Eq. (3.13), is more complicated because we have 
nearest-neighbor interaction (and thus the relevance of 
R,). 



C. Symmetry considerations 



In Sec. HE 1 , after deriving the general Green func- 
tion equation using Q{E), we specialized to rectangular- 
boundary lattices and used lattice reflection symmetries 
to diagonalize the 4x4 matrix Q{E) and obtai ned sc alar 
equations. Here our t-matrix equation Eq. (3.15) re- 
quires us to find the eigenvalues of the t-matrix T . In 
this section, we use particle permutation symmetry and 
pair component exchange symmetry to diagonalize the 
A^'o X A'o t-matrix T{E) for a few special cases. 



1. No = l 

There is only one momentum vector in Qq- Let us write 
Qo = {qi} (this implies that P — qi = qi). Then there 
is only o ne u nperturbed two-particle basis state |qi,qi) 
(see Eq. ( |3.2| )). This must be a boson state, and T{E) is 
a number. We write the resulting scalar equation as, 



E = Eq+Ti{E). 



2. No 



(3.18) 



Here Qq = {qi, q2} with qi +q2 = P. The basis states 
are |qi,q2) and |q2,qi)- The symmetric (boson) combi- 
nation is (|qi, q2) -I- |q2, qi))/'\/2, and the antisymmetric 
(fermion) combination is (|qi,q2) — |q2,qi))/\/2. These 
have to be the eigenvectors of T{E) . And that is to say 
that if we define 



then we have 52 



1 

So, So = /, and 



1 1 
1 -1 



S2T[E)S2- ( ^^'J,^^^ T^XiE) 



(3.19) 



(3.20) 



Here Ti^i{E) and Ti^^i{E) are scalars that correspond to 
boson and fermion symmetries respectively. And our t- 
matrix equation Eq. ( 3.15| ) is reduced to two scalar equa- 
tions, 

E = Eq + T^,i{E), E = Eo + Ti,^iiE), (3.21) 
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for bosons and fermions respectively. Our notation for 
the eigenvalues of T{E) is always to write T with sub- 
scripts that are the coefficients (in order) of the iVo two- 
particle basis vectors. 



3. No =4: 

Here Qo = {qi, q2, qa, q4} with qi -|- qa = qs -I- q4 = 
P. The basis states are |qi,q2), |q2,qi), |q3,q4), and 
|q4,q3). Using particle permutation symmetry, we get 
two states with even symmetries appropriate for bosons, 
which generically would be 

a(|qi,q2> + |q2,qi)) + &(|q3,q4) + |q4,q3)), 

-K|qi,q2> + |q2,qi» +a(|q3,q4> + |q4, q3)C8-22) 
and two odd (fcrmion-type) states, 

a(|qi,q2) - |q2,qi)) + &(|q3,q4> - |q4,q3>), 

-^(|qi,q2) - |q2,qi)) +a(|q3,q4) - |q4,q3{^-23) 

where a and b are arbitrary coefficients to be determined. 

Recall A^o = 4 means the pair (qi,q2) has the same 
total momentum and energy as (q3, q4), which may hap- 
pen for various reasons. When the re ason is the pair 
component exchange symmetry (of Sec. HI A), i.e. qa = 
{qix,q2y) and q4 = {q2x,qiy), then a = h= 1/2, due to a 
hidden symmetry under the permutation 1 3, 2 ^ 4. 
The only effect this permutation has on the momentum 
transfers q^ — q^ is to change the sign of one or both com- 
ponents; but the potential V{v) is is symmetric under re- 
flection through either coordinate axis, hence V{(\i — q^) 
is invariant under the permutation. Since the t-matrix 
depends only on V{cii — q^), it inherits this symmetry. 
Next, if we define 



/fill 




then we have Sa — S\ 



Si = /, 



(3.24) 



and ^4 T(£') ^4 becomes 



diagonal with four eigenvalues of T(E): Ti_i^i^i{E), 

ri,i,_i,_i(i;), ri._i.i._i(^), and ri,_i,_i^i(i;j)'. ^ And 



our t-matrix equation Eq. (3.15) is reduced to 



E = E^+Ti^^,i^i{E), E = Eo+T^,i,-i,-i{E), (3.25) 
for bosons and 

-Ti^-i^i^-i{E), E = Eo + Ti^^i^- 



E — En 



(3.26) 

for fermions. 

The three cases A^o = 1:2, and iVg = 4 with pair 
component exchange symmetry are three special cases in 
which we know the eigenvectors of T and can therefore 
diagonalize T from symmetry considerations easily. 



Different or larger values of Nq are possible when P has 
a special symmetry, e.g. when = Py, Nq = 8 generi- 
cally since Qo includes pairs such as (giy, Oi^, (?2y, q2x)- 
For these general cases, we return to Eq. ( 3.15 ) and diag- 
onalize T numerically. For example, on a i x L lattice, 
the pairs (0,1)(0, —1) and (1,0)(— 1,0) have the same to- 
tal energy and momentum, but this is not due to the pair 
component exchange symmetry. In this case, we numeri- 
cally diagonalize the 4x4 matrix T{E ), and we fi nd that 
in the fermion eigenvectors, Eqs. ( |3.23 ) and (3.23), a ^ b. 



D. Solving for energy 

The example system that we will study here is 10 x 
11 with P = (0,0). The noninteracting and interacting 
energies of the system are in Table I It can be seen that 
all of the energies listed in Table | are of the three cases 



discussed in Sec. IIIC: A^o = 1, Afg = 2, and Afg = 4 due 



to pair component exchange symmetry. 



TABLE I: The 12 low-lying noninteracting and exact two- 
particle energies of the 10 x 11 lattice with total momentum 
P = (0, 0). qi and q2 = P — qi are the momentum vectors. 



qi 


q2 


£:(qi)+£-(q2) 


boson 


fermion 


(0,0) 


(0,0) 


-8.0000000000 


-7.9068150537 


-7.3117803781 


(0,1) 


(0,-1) 


-7.3650141313 


-7.2998922545 


-7.1770594424 


(1,0) 


(-1,0) 


-7.2360679774 


-6.9713379459 


-6.4994071102 


(1,-1) 


(-1,1) 


-6.6010821088 


-6.6010821088 


-6.4700873024 


(1,1) 


(-1,-1) 


-6.6010821088 


-6.0227385416 


-5.5449437453 


(0,2) 


(0,-2) 


-5.6616600520 


-5.4277094111 


-5.1475674826 


(2,0) 


(-2,0) 


-5.2360679774 


-5.0769765528 


-4.8309218202 


(1,2) 


(-1,-2) 


-4.8977280295 


-4.8977280295 


-4.7226011845 


(1,-2) 


(-1,2) 


-4.8977280295 


-4.6571944706 


-4.3808316899 


(2,-1) 


(-2,1) 


-4.6010821088 


-4.6010821088 


-4.1884725717 


(2,1) 


(-2,-1) 


-4.6010821088 


-3.5439149838 


-3.3270813673 


(0,3) 


(0, -3) 


-3.4307406469 


-3.1234645374 


-2.8242092883 



We solve for energy E in the implicit equation, E = 
Ea + T{E), where T{E) represents the eigenvalues of 
T{E), e.g., Ti,_i(£;). We plot f{E) = Eo+T{E) along 
with a line y = E. Their intersections are the desired 
energies E. 



1. No = 1 case 

In Fig. ^, we plot f{E) versus E for the 10 x 11 
lattice with P = (0, 0) and the noninteracting energy 
Eq = -8.0 = E{fd)+E{Q). Here Qo = {(0,0)}, and the 
nonperturbed state is |qi = (0,0),P — qi = (0,0)) which 
can only be a boson state. The energy intersections from 
Fig. I are -7.906, -7.299, -6.971, -6.022, and so on. 
Looking into Table |, we see that these are all boson en- 
ergies. 
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FIG. 3: f{E) ^ Eo + Ti{E) versus E for 10 x 11 lattice 
with P = (0, 0) and Eo = -8.0 (i.e., Eo = £{0) + £{0).) The 
intersections with the line y — E are the exact two-particle 
energies. 



2. No = 2 case 

In Fig. I we plot for Eq = -7.365 and P = (0, 0) 
with Qo = {(0, 1), (0,-1)}. The boson function f{E) = 
Eq + Ti,i{E) is the dotted hne in the top graph, and the 
fermion function f{E) = Eq + Ti^-i{E) is the solid line 
in the bottom graph. 




In Fig. ^, note also that the energy —6.601, which is 
an exact eigenenergy from exact diagonalization, does 
not appear as an intersection in Fig. |3[ This is a special 
energy, being also a noninteracti ng en ergy. Earlier, as 
mentioned at the beginning of Sec. Ill B, we assumed that 
our E ^ £{q) + £{P ~ q) for any q G Q, so this energy is 
excluded from our t-matrix formulation. We will address 
later in Sec. Ill D 3 this kind of exact solutions that are 



also noninteracting energies. 

Note that our equation E = Eq + T{E) is a reformula- 
tion of the Schrodinger equation with certain symmetry 
considerations, and it should be satisfied by all energies 
E with the same symmetry. Building T{E) from Eq and 
P does not automatically give us a unique interacting 
energy E that corresponds to the noninteracting energy 
Eq. However, we can see clearly from Fig. |^, if we per- 
turb the exact solutions by a small amount E E + S, 
then f{E) changes drastically except for the lowest en- 
ergy E = —7.906. That is to say that these other ener- 
gies, for example E = —6.971, are exact solutions of the 
equation f{E) — E, but they are not stable solutions. 
From the plot, only E — —7.906 comes close to being 
stable. 

We can be more precise about this notion of stability. 
If we have an iteration Xn+i = fixn), and x* is a fix point 
(i.e., f{x*) — X*), then the iteration is linear stable at 
X* if and only if |/'(a;*)| < 1. In our plots, we have 
included a line y = E with slope one, which can be used 
as a stability guide. An intersection (fix point) is linearly 
stable when the function f{E) at the intersection is not 
as steep as the straight line. 




FIG. 4: f(E) versus E for 10 x 11 lattice with P = (0, 0) and 
Eo = -7.365 (i.e., Eo = £{0, 1) + £{0, -1).) The top graph 
(dotted line) is for boson f{E) — Eo + Ti,i(_B), and the top 
graph (sohd line) for fermion f{E) = Eo +Ti.-i{E). The 
fermion curve is essentially flat near E = Eo- 

The intersections closest to Eq — —7.365 are —7.299, 
the first excited boson energy (see Table ^, and —7.311, 
the lowest fermion energy. Note that the curve on which 
the fermion intersection (—7.311) lies is very flat. In 
other words for this fermion energy E k Eq + T(Eq), i.e., 
the first iteration using the noninteracting energy gives 
an energy very close to the exact value. More precisely, 
we find with Eq = -7.365014, /{Eq) = Eq + T{Eq) = 
-7.310584, which is vep,|Glo»p to £: = -7.31178. Many 
t-matrix calculations ,E2lESEal23 use the first iteration E w 
Eq + T(Eq) as an approximation to the exact energy, and 
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we see in this case this approximation is very g ood. (We 
wih come back to this point later in Sec. Ill E .) 



3. No = A case 
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FIG. 5: Boson f{E) versus E for 10 x 11 lattice with P = 
(0,0) and Eo = -6.601. The dotted line is for ri,i,i,i and the 
horizontal dot-dashed line for Ti,i,_i,_i (which corresponds 
to a noninteracting state, see text at the end of this section). 
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-9.0 ^ 
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FIG. 6: Fermion f{E) versus E for 10 x 11 lattice with 
P = (0,0) and Eo = -6.601. The solid line is for Ti,_i,i _i 
and the long-dashed line for Ti^-i^-i,!. Note that closely 
spaced fermion energy pairs are separated by symmetry. 



In Figs. I and |, we plot f{E) for E^^ = + 
£(-1,1) = £(1,1) +£(-1,-1) = -6.601 and P = (0,0). 



For this A*o = 4 case we have two boson functions, plot- 
ted in Fig. |, f{E) = Eo + Ti^i^i^i{E) (dotted line) and 
f{E) — £'o+ri4^_i^_i(i?) (dot-dashed line) , and we have 
two fermion functions, plotted in Fig. ^ f{E) — Eo + 
Ti^-iX-i{E) (solid line) and f{E) = Eq + Ti^^i^^i^i{E) 
(dashed line) . The fermion intersections closest to Eq are 
—6.499 and —6.470. Here again the two fermion curves 
are very flat. The two boson intersections closest to Eq 
are —6.022 and —6.601. Note that the latter is also a 
noninteracting energy, and it is the intersection of the 
horizontal line E = Eq with y = E. 

One interesting observation of the fermion plot in 
Fig. H is that pairs of closely spaced energies (for example 
— 7.311 and —7.177) lie on different symmetry curves. We 
know that if we have a square lattice (for example 10 x 10) 
then the noninteracting fermion energies come in pairs. 
Here, we have chosen a 10 x 11 lattice that is close to 
a square but does not have exact degeneracies. We see 
that the resulting closely spaced pairs are separated by 
symmetry considerations. 

Another interesting observation from Fig. ^ for bosons 
is that we have a horizontal line that corresponds to 
Ti,i__i^_i(i?) = 0. For this case the noninteracting en- 
ergy is an exact energy. That is to s ay, (1,1,-1,-1) is a 
null vector of T{E) (see Sec. IIIC3), or the eigenstate. 



|qi,q2) 



|q2,qi; 



|q3,q4) 



|q4,q3; 



(3.27) 



with qa = {qix,q2y) and q4 = {q2x,qiy) is an ex- 
act eigenstate of the Hamiltonian. This can be shown 
easily using the Schrodinger equation Eq. ( p. 15 ). We 
have 5r(qi) = ^(qa) = 1, ^(qs) = 5(q4) = -1, and 
g{q) = for all other q, and we can easily show 
V{q ~ qi) + V{q - qa) - V{<1 - qg) - ^(q - q4) = 
(because V^(k) can be separated into a sum of two terms 
that involve the x and y components separately). 

Transforming to the real space, without worrying 
about normalization, we can have 



g(r) = Y^e-^^'^gid) 



where we have used the fact mentioned above that (?(q) is 
not zero for only four q's which are related by pair com- 
ponent exchange symmetry. It is clear from Eq. ( 3.2S| ) 
that 5(0,7/) = = 5(2;, 0), which means that the wave 
function in relative position space is "d-wave" like, hav- 
ing nodes along x and y axes (thus happens to have nodes 
at every relative position where the potential would be 
nonzero). 



E. Fermion: noninteracting to interacting 

In this section we use the t-matrix techniques devel- 
oped in the preceding sections of this section to study 
the relationship between the noninteracting energies and 
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the interacting energies. We start with the table of ener- 
gies in Table | for the 10 x 11 lattice with P = (0, 0). We 
have asked in the introduction to this section whether 
we can go from the noninteracting to the interacting 
energies and now we know that we have an equation 
E = Eo + T{E) where T{E) is the symmetry reduced 
scalar t-matrix function. From our graphs (Fig. ^ and 
Fig. ^) we have commented that for fermions the curve 
of T{E) around Eq is quite flat (which is not the case for 
bosons) . And we mentioned that this implies that the ap- 
proximation E « Eq+T{Eq) is close to the exact energy. 
Now in this section, we study the t-matrix approach for 
a specific system. We will denote Ei = Eq + T{Eo), the 
first iteraction result, and £'„+i — Eq + T{En), the nth 
iteration result. 

In Table || we show the t-matrix calculation for the 
10 X 11 lattice. We show for the lowest few states the non- 
interacting energy Eq, the first t-matrix iteration Ei, the 
fifth t-matrix iteration E's, and the exact energy £^cxact- 
In Fig. these energy levels are plotted graphically. From 
the table, it is clear that the first t-matrix iteration re- 
sult El is quite close to the exact energy, and the fifth 
iteration result E^ gives a value that is practically indis- 
tinguishable from the exact value. 



EO 



E1 



-5.0 



-6.0 



-7.0 



E5 



Exact 



i.O L 



-5.0 



-6.0 



-7.0 



-8.0 



ground state energy of a few particles on a large lat- 
tice by summing up the energy of each pair of particles. 
This section contains a much more detailed study of the 
few-fermion problem: we will consider first the fermion 
shell effect and then we will study the interaction correc- 
tion to energy (ground state and excited states) for a few 
fermions (three, four, and five) using the t-matrix. 
Our results - summarized in Sec. IV F - confirm that. 



in the dilute limit, almost all of the interaction correction 
is accounted for by the two-body terms of the t-matrix 
approximation, Eq. (4.1). But (recall Eq. (1.3)) that is 
a hallmark of a Fermi liquid picture; i.e., our numerical 
results suggest its validity at low densities. This is a non- 
trivial result, in that firstly, the validity of Fermi liquid 
theory in a finite-system context has rarely been con- 
sidered. Standard t-matrix theory depends on a Fermi 
surface which (at T = 0) is completely sharp in momen- 
tum space, and every pair's t-matrix excludes scattering 
into the same set of occupied states. In a finite system, 
however, the allowed q vectors fall on a discrete grid, 
and since the total number of particles is finite, the t- 
matrices of different pairs see a somewhat different set 
of excluded states (since they do not exclude themselves, 
and one particle is a non-negligible fraction of the total). 

Secondly, and more essentially, the analytic justifica- 
tions of Fermi liquid theory exist only in the cases of spin- 
full fermions (in a continuum). That case is dominated 
by s-wave scattering, so that the t-matrix approaches a 
constant in the limit of small momenta (and hence in the 
dilute limit). Our spinless case is rather different, as will 
be elaborated in Sec. because the t-matrix is dom- 
inated by the p-wave channel, which vanishes at small 
momenta. Thus the q dependence is crucial in our case, 
and the numerical agreement is less trivial than it would 
be for s-wave scattering. 

In this s ection, after an exhibition of the shell effect 
(Sec. IV A), we present a general recipe for the multi- 
fermion t-matrix calculation. This is developed by the 
simplest cases, chosen to clarify when degeneracies do or 
do not arise. 



Fermion shell effect 



FIG. 7; Two-fermion energy levels for the 10 x 11 lattice with 
P — (0,0). From left to right, the lowest few noninteracting 
energies Eq, first t-matrix iteration Ei, fifth t-matrix itera- 
tion Es, and the exact are plotted. Note that 
the third noninteracting energy from the bottom is doubly 
degenerate (see Table O). 



IV. A FEW FERMIONS: SHELL EFFECT AND 
T-MATRIX 

In Sec. ||, we used lattice Green function to study the 
problem of two particles (bosons and fermions), and at 
the end of that section, in Sec. [ID, we obtained the 



At zero temperature, the ground state of noninteract- 
ing fermions is formed by filling the one-particle states 
one by one from the lowest to higher energies. For our 
model of spinless fermions on a square lattice, we have the 
two ingredients for the shell effect: fermionic exclusion 
and degeneracies of one-particle states due to the form 
of our energy function and lattice symmetry. Shell effects 
have been noted previously in interacting models ;c3 our 
code, permitting non-rectangular boundary conditions, 
allows us to see even more cases of them 

In Fig. ^ we show the exact and for comparison the 
noninteracting ground state energies for the 5 x 8 and 7x7 
lattices for up to seven particles. The energy increment 
curve E{M) ~ E{M — 1) is plotted and shows clearly the 
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TABLE II: Fermion energies for 10 x 11 lattice with P = (0,0). -Eo = f (qi) + £{<l2) is the noninteracting energy. En — 
Eo + T{En-i) where T{E) is the symmetry reduced t-matrix. Here only fermion energies (from Ti,_i or Ti.-i,!,-! and 
ri,_i,_i,i) are included. 



qi q2 -Eo El E5 -Ecxact 

(0,1) (0,-1) -7.365014 -7.310598893 -7.311780378 -7.311780378 

(1.0) (-1,0) -7.236067 -7.17521279 -7.17705944 -7.177059442 
(1,-1) (-1,1) -6.601082 -6.493807907 -6.49940706 -6.49940711 

(1.1) (-1,-1) -6.601082 -6.460962404 -6.470087137 -6.470087302 
(0,2) (0,-2) -5.661660 -5.532751985 -5.54494225 -5.544943745 
(2,0) (-2,0) -5.236067 -5.134290466 -5.147558003 -5.147567483 



shell effect. ison, we show the boson energy plot for the 5x8 lattice in 

Fig. ^. Because bosons can all be at the zero-momentum 
state, where energy is —4, the total noninteracting energy 
I • • • I is —AM. The exact energy curve shows smooth changes 
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FIG. 9: Boson noninteracting and exact groundstate energies 
for the 5x8 lattice with one to seven particles. Because bosons 
can all be at the zero-momentum state, where the energy is 
—4, the noninteracting energy is —4M. The exact energy 
shows smooth changes when M increases. There is no shell 
effect. 
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FIG. 8: Shell effect for 5 x 8 and 7x7 lattices. Exact, 
interacting groundstate energies are compared with noninter- 
acting energies for up to seven particles. Energy increment 
E{M) - E{M - 1) is shown. 



The filled shells for the 5x8 lattice are at M = 3 (with 
momentum vectors (0, 0)(0,±1) occupied) and M = 5 
(with (0, 0)(0,±1)(±1,0) occupied). On the other hand, 
M = 3 is not a filled shell of the 7x7 lattice. For compar- 



B. General multi- fermion theory 

The key notion for generalizing our two-ferniion ap- 
proach to A/ fermions is that the set Qq now consists of 
every M-tuple a of wavevectors that gives the same total 
momentum and noninteracting energy. This defines a re- 
duced Hilbert space, with the corresponding basis states 
\4>a)- We can construct an approximate, effective Hamil- 
tonian Hq + Htm acting within Qo-space, where Htm is 
a sum of pairwise t-matrix terms, each of which changes 
just two fermion occupancies: 



Ht, 



-£1/3 



(4.1) 
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The notation X^l/? means the sum only includes the pair 
{a, (3) when \cj)a) differs from jt/i/s) by a change of two 
fermions. 



States 



Thus, each term in Eq. (41) is associated with a par- 
ticular fermion pair (q^ , C[j ) . Each such pairwise t-matrix 
can be viewed as a sum of all possible repeated scatterings 
of those two fermions through intermediate states, ex- 
cept that intermediate states which are already included 
in Qq are excluded. (The most important omission of 
this approximation would be the processes in which three 
or more fermions are permutated before the system re- 
turns to the Qo Hilbert space.) Each term is a two- 
fermi on t-m atr ix cal culated according to the recipe of 
Sees. [II B and [II D| - thus each term has its own two- 
fermion wavevecto r se t Qo''"' aud complementary set Q^'^ , 
as defined in Eq. (p]). The only change in the recipe is 
to augment the set Q*'^ of wavevectors forbidden in the 
intermediate scatterings of the two fermions, since they 
cannot scatter into states already occupied by the other 
fermions in states a and (3. (See Eq. (4.2) for an exam- 
pie.) 

The t-matrix treatment is a form of perturbation ex- 
pansion, for which the small parameter is obviously not 
V (which is large) but instead 1/i^, as is evident from 
Eq. (2.44). That is, as the lattice size is increased (with 
a fixed set of fermions), the approximation captures a 
larger and larger fraction of the difference i?oxact — Eq. 



C. A three-fermion t-matrix calculation 

We first compute the energy of three fermions (Af = 3) 
for the 8x9 lattice with P = 0. For this example calcu- 
lation, we have chosen ^ Ly to reduce the number of 
degeneracies in the noninteracting spectrum. In Fig. ^ 
we show the lowest five noninteracting levels and the cor- 
responding states in momentum space. 

Let us consider the lowest noninteracting state in the 
8 X 9, P = (0, 0), and M — 3 system, with three momen- 
tum vectors: qi = (0,1), q2 = (0,0), and qa = (0,-1) 
(see Fig. |l^). And let us first consider the interaction 
of the pair qi and q2. The noninteracting energy of the 
pair is E^"^ = S (qi) + £ {c[2) = -7.682507 and the total 
momentum is P12 = qi + q2 = (0, 1). As usual, we use 
and P12 to form the set Qo^^ (Eq. (3.1)). Here there 



Eq . , _ , ,,, 

are no other degenerate vectors so Qo^^ — {(0, 0), (0, 1)}. 
The three-particle problem is different from the two- 
particle case in the choice of Q^^, the set of momentum 
vectors that the two particles can scatter into. Due to the 
presence of the third particle and Pauli exclusion, the two 
particles at qi = (0,1) and q2 = (0,0) cannot be scat- 
tered into the momentum vector qs = (0, —1), so we must 
exclude qs from Q^^. Furthermore, even though there is 
no particle at P12 - qs = (0, 1) - (0, -1) = (0, 2), this 
momentum cannot be scattered into, because otherwise 
the other particle would be scattered into the occupied 
qa. That is to say, the momentum vectors that can be 



Energy 
-11.06417777 

-10.82842712 
-9.892604897 

-8.694592710 
-8.239901252 



FIG. 10: Lowest five noninteracting energy levels for the 8x9 
lattice with M = 3 fermions and total momentum P = (0, 0). 
States in momentum space are drawn. 




scattered into are 



= {q|q 7^ qi,q2,q3,Pi2 - qa} 



(4.2) 



This exclusion is shown graphically in Fig. |Tl|. 

The t-matrix formalism can then be applied using Qo^^ 
and Q^^ to compute energy correction T^'^{Ei2) for the 
interaction of the qi and q2 pair. Here T^^JE) here 
is the "fermion" function Ti_i{E) (Eq. (|3.2l|) ), cor- 
responding to the antisymmetric eigenvector of the t- 
matrix 'T{E); the tilde denotes the modification due to 
exclusion of the set Q^^. When the t-matrix contributes 
a small correction, it is accurate to use the bare values, 
Efj = £{(ii) + £{qj), and this approximation was used 
for all tables and figures in this section. 

The total energy within this approximation is a sum 
of the t-matrix corrections for all possible pairs in the 
system, which are (q2,q3), and (qi,q3) in the present 
case: 



Etm = £{cii)+£{q2)+£iq3)+T''iEi2)+T'%Ei3)+T'\E23) 

(4.3) 

This is a special case of the effective Hamiltonian 
Eq. (4.1), which reduces to a 1 x 1 matrix in the non- 



degenerate case. (That is, whenever the set Qq of multi- 
fermion occupations has just one member.) The momen- 
tum space exclusions due to the presence of other parti- 
cles are depicted in Fig. |l^, and the numerical values of 
this calculation are given in Table III. 

A more accurate approximation is to enforce a self- 
consistency, 



Ej. 



E^^+T={E,j) 



(4.4) 



where as defined above i?° = f(qi) +£{cij). It should 
be cautioned that the physical justification is imperfect: 
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if we visualize this approximation via a path integral or 
a Feynman diagram, the self-consistent formula would 
mean that other pairs may be scattering simultaneously 
with pair (ij), yet we did not take into account that the 
other pairs' fluctuations would modify the set of sites 
Q^^ accessible to this pair. In any case, analogous t o the 
two-particle t-matrix (Sec. Ill), we could solve Eq. (4.4) 
iteratively setting E^+^ = Ef^ +f'3{E^j), until succes- 
sive iterates agree within a tolerance that we chose to be 
10~^^, which happened after some tens iterations. 




FIG. 11: Momentum space exclusions in t-matrix M = 3 
calculation for the state (0,0)(0,1)(0,-1). The crosses indicate 
exclusions when calculating pair energy for (0,0) (0,1) (left fig- 
ure), (0,0) (0,-1) (middle), and (0,1) (0,-1) (right) respectively. 



TABLE IV: Lowest 15 noninteracting, exact, and t-matrix 
energies for 8x9 lattice with M — 3 and P — (0, 0). 
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TABLE III: T-matrix calculation for the 8x9 lattice with 
M — 3 noninteracting particles qi — (0,0), q2 ~ (0, 1), and 
qa — (0, —1). The total noninteracting energy is Eq — f (qi) + 
£(q2) + £{(13) and the total t-matrix correction is T = T^^ -|- 
j,i3 j_ y23 rpj^g energy calculated using the t-matrix is then 
Etm = Eq + T and the exact energy from diagonalization is 
inexact. E'g = £{(ii) + £{<lj), is the noninteracting energy of 
the pair. 





P 




rpij 


(0,0)(0,1) 


(0,1) 


-7.532088886 


0.041949215 
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(0,-1) 


-7.532088886 
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Column sum 
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= 0.202583012 


Noninteracting 


total 


Eo = 


-11.064177772 


T-matrix total 




-Etm = 


-10.861594761 


Exact total 




-Ecxact — 


-10.871031687 



Using the same procedure, we can also calculate the 
t-matrix energies for the nondegenerate excited states 
of the M = 3 system in Fig. |ig: the (-1, 0)(0, 0)(1, 0) 
and (0, 2)(0, 0)(0, —2) states. The results are shown in 
Table IV. Fig. |lj shows graphically the noninteracting 
energy levels, the t-matrix energies for the three nonde- 
generate states, and the exact energies from diagonaliza- 
tion, and the arrows link the noninteracting energies i?o 
with the t-matrix results ii^tm = £-0 + T. The agreement 
between £^tm and -Eoxact is good. 



D. A five-fermion t-matrix calculation 

We now consider briefly a M = 5 calculation, again 
for the 8x9 lattice. The noninteracting ground state is 
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FIG. 12: Noninteracting, t-matrix, and exact energies of 
the three-particle fermion system on the 8x9 lattice with 
P = (0,0). The bracketed numbers refer to the degeneracies 
of the level (see Fig. px| ). The arrows associate the noninter- 
acting states with the t-matrix results. We have worked on 
nondegenerate noninteracting states so far. 



unique, with momentum vectors qi = (0, 0), q2 = (0, 1), 
qs = (0,-1), q4 = (1,0)_, andqs = (-1,0). In Fig. |l| we 
show the excluded set Q^'^ of the t-matrix computation 
for the pair (q2, q4)- The momentum vectors (qi, qs, qs) 
filled with other fermions are excluded, of course; three 
more wavevectors are excluded since the other member 
of the pair would have to occupy one of qi , qa , or qs , due 



17 



to conservation of the total momentum P = (1, 1). The 
t-matrix results for all 10 pairs are presented in Table 0. 

One might think that the pair, q2 = (0, 1) and q4 = 
(1,0), exhibits pair-excha nge sym metry w ith (0,0 )(1,1), 



so that A'o = 4 as in Sec. Ill C 3 and Sec 



HID 3. How- 



ever, since (0,0) is occupied, the (0,1)(1,0) pair cannot 
scatter into the (0,0)(1,1) pair: hence (0,1)(1,0) is a 
generic pair with A^q'* = 2. In general, if a pair is ever 
free to scatter into a degenerate pair state with a different 
occupation, that must be part of a many-particle state 
degenerate with the original one. Thus, the complicated 
t-matrix pairs with Nq > 2 can arise in a many-fermion 
calculation only when the noninteracting many-fermion 
states are themselves degenerate. 




FIG. 13: Momentum space exclusions in t-matrix AI = 5 
calculation for momentum vectors q2 — (0, 1) and q4 — (1,0) 
(dots without crosses). These two fermions are excluded 
from scattering into momenta from the set Q^'* (marked by 
crosses). The ground state is shown, with occupied momenta 
(0,0), (0,1), (0,-1), (1,0), and (-1,0) (solid dots). 



TABLE V: T-matrix 
five particles (0,0), (0 
exclusions in Q for the 



calculation for the 8x9 lattice with 
,1), (0,-1), (1,0), and (-1,0). The 
pair (0,1) (1,0) are depicted in Fig. hA. 
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-7.532088886 


0.045184994 


(0,0) 


0,-1) 


-7.532088886 


0.045184994 


(0,0) 


1,0) 


-7.414213562 


0.056898969 


(0,0) 


-1,0) 


-7.414213562 


0.056898969 


(0,1) 


0,-1) 


-7.064177772 


0.118684581 


(0,1) 


1,0) 


-6.946302449 


0.081095408 


(0,1) 


-1,0) 


-6.946302449 


0.081095408 


(0,-1 


) (1,0) 


-6.946302449 


0.081095408 


(0,-1 


) (-1,0) 


-6.946302449 


0.081095408 


(1,0) 


-1,0) 


-6.828427125 


0.131405343 




Eo = 


-17.892604897 


f = 0.778639481 




-E'cxact — 


-17.145715214 £tm 


= -17.113965417 



E. Degenerate states 

In the ground state examples considered up to now 
(Sees. IV C and IV D), the noninteracting states were 
Let us now study a degenerate state 

3 



all nondegenerate. 

in the third lowest level (six-fold degenerate) of M 



fermions on the 8x9 lattice : q2 = (0, 1), qa = (1,0), 
q4 — (—1,-1). (See Fig. |l0|, row 3.) In this state, the 
pair [q2 = (0,1), qa = (1,0)] has the same total energy 
and momentum as the pair [qi = (0,0), qs — (1, 1)], on 
acco unt o f the pair component exchange symmetry (see 
Sec. [II A ); consequently [q2,q3l can be scattered into 



[qi,q5] contrary to the previous example in Sec. IV E . 
Indeed, each of the six basis states in row 3 of Fig. |10| 
is connected to the next one by a two-body component 
exchange symmetry. 

Following th e two-f ermion calculation with A^o = 4 
pairs (see Sec. IIIC 3), the degenerate pairs q2q3 and 



[qi,q5] mustbe handled in the same set Qo^^- The re- 



suhs Eqs. ( |3^ ), (|^, and ( ^^ imply 



T"''(£^23)|q2q3) 
1 



1 



(Ti,_i,i,_i-t-ri,_i,_i,i)|q2q3) 



+ 2 ( Ti,_i,i,-i-Ti,_i,_i,i)|qiq5)(4.5) 

Here Ti,_i4,_i and Ti^-i^-i^i depend implicitly on P = 
(0, 0), on the momenta, and on the energy -^2,3, as well as 
on Q^^ which depends on the occupation (q4) of the third 
fermion. In this notation, each T'' -^ acting on any state 
produces two terms as in Eq. (4.5). The total t-matrix 
correction Hamiltonian is J^ij T''-* , summed over all 18 
possible pairs appearing in the degenerate noninteracting 
states. When we apply this to each state in the third 
row of Fig. 10, we finally obtain a 6 x 6 matrix mixing 
these states. Diagonalization of this matrix would give 
the correct t-matrix corrections (and eigenstates) for this 
"multiplet" of six states. We have not carried out such a 
calculation. 

It is amusing to briefly consider the states in row 5 
of Fig. |l^, a different sixfold degenerate set. Unlike the 
row 3 case, these states separate into two subsets of three 
states, of which one subset has {qy} — —2,-1-1,-1-1 and 
the other subset has the opposite qy components. Scat- 
terings cannot mix these subsets, so the 6x6 matrix 
breaks up into two identical 3x3 blocks. Hence the 
t-matrix energies from row 5 consist of three twofold de- 
generate levels. By comparison, the exact interacting en- 
ergies derived from these noninteracting states come in 
three nearly degenerate pairs, such that the intra-pair 
splitting is much smaller than the (already small) split- 
ting due to the t-matrix. 



F. Errors of the t-matrix 

How good are the t-matrix results? From our example 
calculations on the 8x9 lattice, in Tables HI, IV, and 0, 
we see that Etm and -Ecxact are close. 

In Fig. |lj we plot the noninteracting, t-matrix, and ex- 
act energies for M = 3, P = (0, 0) ground state on a se- 
ries of near square lattices Lx (L+l). The noninteracting 
ground state momentum vectors are (0, 0)(0, 1)(0, —1) for 
this series of lattices. We do not plot for L > 12, be- 
cause, as can be seen in the bottom graph, the t-matrix 
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energy iJtm approaches the exact energy i^oxact rapidly. 
To see more clearly the error of the t-matrix result, we 
plot also Etm — £^cxact, which decays very fast as the size 
of the lattice increases - very roughly as the power. 
Even at L = 6, i.e. at a density n « 0.05, the t-matrix 
approximation captures 95% of the interaction energy 
Ecxact — Eq. These figures are based on using the bare 
energies in T^\Eij ) in Eq. ( |4.3| ). If we carried out the 



self-consistent calculation described in Sec. IV C), the er- 
ror i?tm — Ecxact would bc Smaller by a factor of roughly 
2.5. 




FIG. 14: Noninteracting, t-matrix, and exact energies for 
M = 3, P = (0,0) ground state (0,0)(0,1)(0, -1) on a series of 
Lx (L + l) lattices as a function of L (top graph). i?tm — £^exact 
versus L (bottom graph). 



V. THE DILUTE LIMIT: ENERGY CURVES 

In this section, we are interested in the functional form 
of the energy as a function of particle density for both 
bosons and fermions in the dilute limit. In the three- 
dimensional case, the problem of dilute quantum gases 
with strong, repulsive, short-range interactions was first 
addressed ip_.the language of diagrammatic field theory 
by Galitskiio for fermions and BeliaevcZl for bosons. At 
that time, the ground state energy as an expansion in 
the particle density was also obtained for hard-sphepi 
fermion and boson gases by Yang and coUaboratorsea 
using a pseudopotential method. The field theoretical 
methods werCrlater adapted to two dimensions in pMlic- 
ular by SchickEll for hard-disk bosons and by Bloomc3 for 
hard-disk fermions. Some other relevant analytic papers 



using a t-matrix approach for t. 
discussed in Sec. [ B 



Hubbard i 
t^mamoriCJ and Mattis 
for d = 2. 



pdel were 
I in d = 3 

and Rudin and Mattia 

For both hard-disk fermions and bosons in two dimen- 
sions, the leading-order correction to the noninteracting 
energy is found to be in the form of n/lnn, where n 
is particle density. Expansions with second-order coef- 
ficients different from the results of Schick and Bloom 
were found in Refs. ^ and |5^ for the boson case and in 
Refs. |4| and ^ for the fermion case. There is no con- 
sensus at this time on the correct second-order coefficient 
for both the boson and fermion problems. 

Recently, Ref. ^ has proved rigorously the leading- 
order expansion ptthe two-dimensional dilute boson gas 
found by Schick.EiJ Numerically, the dilute boson prob- 
lem on a two-dimensional lattice has been studied us- 
ing quantum Monte Carlo in Refs. ^ and ^ and they 
obtain good fit with Schick's result. As we mentioned 
in Sec. IB, more recently, because of a question re- 



garding the validity of the Ferpaj liquid theory in two 
dimensions- Bloom's calculationEj has received renewed 



attention! 



numerical studies. 



but this result has not been checked by 



Dilute bosons 



For two-dimensional hard disk bosons, the energy per 
particle E/M at the low-density limit from diagram- 
matic calpiilations is obtained (in the spirit of Ref. |4^) 
by SchickEl 



E 2TTh-' 



M m |ln(ria^)| 



l + O 



1 



ln(na^ 



(5.1) 



where n = M/N is particle density, m the mass of 
the boson, and a the two-dimensional scattering length. 
As mentioned above, the coefficient of the second-order 
term, has not been settled. 

This hard-disk calculation was carried out using the ki- 
netic energy h^k"^ /2m. In our l attice model, our hopping 



energy dispersion is (Eq. ( ^.11 )) 



£(k) = —2t{cos{kx) + cos{ky)) 



-At + tk^ 



(5.2) 



where we have Taylor-expanded the dispersion function 
near k = because in the dilute limit, at the ground 
state, the particles occupy momentum vectors close to 
zero. Therefore if we use t = h = 1 and the effective 
mass m* such that we have the form U^k'^ /2m*, then 
TO* = 1/2 for our system. So for our model, Schick's 
expansion Ec^. (5.1) should become, 



E 
M 



= -4- 



47rn 



I ln(na*^)| 



l + O 



1 



ln(na*^) 



(5.3) 



where we have used a* to denote the scattering length 
in our lattice system. There is no straightforward corre- 
spondence between Schick's scattering length a in the 
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continuum and our a* on the lattice. With infinite 
nearest-neighbor repulsion, the closest distance that our 
particles can come to is \/2. We expect roughly 1 < a* < 
\/2, and will determine a more precise value from curve 
fitting. 

In Fig. ^ we show the boson energy per particle 
(E/M) versus particle per site (M/N) curve for ten lat- 
tices, ranging from 25 sites to 42 sites, with three or more 
particles {M > 3). The data from all these lattices col- 
lapse onto one curve, especially in the low-density limit. 
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FIG. 15: Boson energy per particle E/M versus particle 
density M/N data for ten lattices and M > 3. Data from 
different lattices collapse onto one curve. The solid line corre- 
sponds to the fitting function —4:+4Tm{A+B/\ ln(na*^)|) with 
a* = 1.36, A = -0.016, and B = 0.959, which is Eq. (0) 
with parameters from Table 



Eq. (5.3), Schick's result apphed to our model, suggests 
the following leading order fitting form for E/M versus 
n at the low-density limit. 



E/M - 



Airn 



= A 



B 



ln(na*^)| 



(5.4) 



That is to say, if we plot {E/M + 4)/(47rn) versus 
1/1 ln(na*^)|, then, if Schick is correct, we should get a 
straight line, with intercept A = and slope B — 1, with 
one adjustable parameter a*. 

In Fig. m, we plot {E/M + 4)/(47rn) versus 
1/1 ln(na*^)| for the low-density limit {n < 0.15) for three 
choices of a* = 1.0, 1.36, \/2. The data points appear to 
lie on straight lines. For a* — 1.36 the fitted intercept is 



A = -0.016 and the slope B = 0.959. In Table [VJ we 
show the fitted slope and intercept for a number of a* 
choices. The slope is zero close to a* = 1.34 and the in- 
tercept is zero close to a* = 1.39. Our data thus suggest 
a* = 1.36 ±0.03. 

5, the solid line is the hmction —4 + 4:7m{A + 
^1) using a* = 1.36, A = -0.016, and B = 
0.959, and we obtain a good fit up to n = 0.15. 



In Fig. 
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FIG. 16: (£/M 4- 4)/(47rn) versus l/|ln(na*^)| plot to 
check Sch ick's formula for two-dimensional dilute bosons 
(Eq. (p!3|)). The data points are for M > 3 and n < 0.15 
from those in Fig. ^ for lattices from 5 x 5 to 6 x 7. For the 
three a* values, the a* — 1.36 choice gives A — —0.016 ~ 
and B = 0.959 « 1. 



TABLE VI: Intercept A and slope B in linear fitting (E/M + 
4)/(47rn) versus l/|ln(na*)| for bosons, using Eq. (5.4). The 
slope is one close to a* = 1.33 and the intercept is zero close 
to a* — 1.39. So we get a* = 1.36 ±0.03. The fitting for three 
choices of a* is plotted in Fig. |l^. 



a* B A 


a* B A 


a* B A 


1.00 1.855 -0.251 
1.10 1.547 -0.178 
1.20 1.289 -0.112 

1.30 1.072 -0.050 

1.31 1.053 -0.044 


1.32 1.033 -0.039 

1.33 1.014 -0.033 

1.34 0.995 -0.027 

1.35 0.977 -0.021 

1.36 0.959 -0.016 


1.37 0.941 -0.0099 

1.38 0.923 -0.0043 

1.39 0.906 0.0013 

1.40 0.889 0.0069 
1.414 0.865 0.015 



For bosons, quantum Monte Carlo can be used to 
obtain zero temperature energies for reasonably large 
systems. For a dilute boson gas on a square lattice 
with on-site hardcore but not nearest-neighbor interac- 
tion, Ref. |5^ has fitted the first term of Schick's formula 
Eq. (|]l), and Ref. ^ has used higher-order terms and 
included the fitting of the chemical potential also. The 
agreement is good in both studies. 



B. Dilute fermions 

For fermions, it customary to write the energy per par- 
ticle expansion in terms of the Fermi wavevector kp. For 
two-dimensional dilute hard disk fermions with a general 
spin s, the energy per particle from diagrammatic calcu- 
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lations, is obtained by Bloo: 



mi 




(see Ref. ^for the spin-1/2 calculation and Ref. 54 for 
general s). 

Eq. means that for our spinless fermions {s — 0), 



the leading order correction to the noninteracting en- 
ergy in Eq. (5.5) is zero, which is due to the fact that 
Eq. (5.5) is derived for s-wave scattering. In our model, 
without spin, only antisymmetric spatial wavefunctions 
are allowed for fermions, and therefore the leading-order 
correction to the noninteracting energy should be from 
p-wave scattering. Ref. ^ contains a formula for p-wave 
scattering in three dimensions where the leading-order 
correction to E — Eq is proportional to {kpa)^ while the 
s-wave correction is proportional to kpa. We are not 
aware of ^ two-dimensional p-wave calculation in the 
literaturej£j and we have not worked out this p-wave 
problem in two dimensions. We expect that the p-wave 
contribution to energy should be considerably smaller 
than that from the s-wave term. In Sec. II D, we have 
considered the case of a few fermions on a large L x L 
lattice, and in Fig. |l| we have studied the interaction cor- 
rection to the noninteracting energy AE. It was shown 
there that AE for our spinless fermions is much smaller 

than that for bosons. 

Using kp 



47rri/(2s -f- 1), we can rewrite Eq. (5.5) as 



E 
M 




(2s + l)m 



In t his form, it is revealed that the second term of 
Eq. ( ^.6[ ) is id enti cal to the first term of the boson ex- 
pression Eq. (5.1), apart from the replacement n — > 
2sn/(2s + 1). In other words, the dominant interac- 
tion term for spinfuU fermions is identical to the ferm 
for bosons, provided we replace n by the density of all 
spin species but one, i.e. of the spin species which can 
s-wave scatter off a given test particle. 



VI. CONCLUSION 

We have studied a two-dimensional model of strongly- 
interacting fermions and bosons. This model is the sim- 
plest model of correlated electrons. It is very difficult to 
study two-dimensional quantum models with short-range 
kinetic and potential terms and strong interaction. There 
are very few reliable analytical methods, and many nu- 
merical methods are not satisfactory. With our simplified 
model of spinless fermions and infinite nearest-neighbor 
repulsion, we can use exact diagonalization to study sys- 
tems much larger (in lattice size) than that can be done 
with the Hubbard model. One of our goals is to publicize 
this model in the strongly-correlated electron community. 



In this paper, we made a systematic study of the dilute 
limit of our model, using a number of analytical tech- 
niques that so far have been scattered in the literature. 
We studied the two-particle problem using lattice Green 
functions, and we demonstrated the use the lattice sym- 
metry and Green function recursion relations to simplify 
the complications brought by nearest-neighbor interac- 
tions. We derived in detail the two-particle t-matrix 
for both bosons and fermions, and we showed the dif- 
ference between the boson and fermion cases and that 
for fermions the first t-matrix iteration is often a good 
approximation. We applied the two-fermion t-matrix to 
the problem of a few fermions, with modifications due to 
Pauli exclusion, and showed that the t-matrix approxi- 
mation is good for even small lattices. 

It is somewhat puzzling that with the essential role the 
t-matrix plays in almost every calculation in the dilute 
limit with strong interactions, no systematic study of the 
t-matrix for a lattice model has been made, as far as we 
know. We believe that our work on the two-particle t- 
matrix and the few-fermion t-matrix is first such study. 
Some approximations that are routinely made in t-matrix 
calculations are graphically presented, especially the use 
of first t-matrix iteration in calculating fermion energy. 
And we demonstrate the qualitative difference between 
the boson and fermion t-matrices. We believe that this 
study is a solid step in understanding dilute fermions in 
two dimensions, and is of close relevance to the 2D Fermi 
liquid question. 

The dilute boson and fermion energy per particle 
curves were studied in Sec. The boson curve was 
fitted nicely with a previous diagrammatic calculation, 
and our work on dikijte bosons complements quantum 
Monte Carlo results S For the fermion problem in our 
model, the leading order contribution to energy is from 
p-wave scattering; therefore, the_series of pCESults based 
on s-wave calciilations by Bloomj£3 Bruch,EJ and Engel- 
brecht, et alEB are not directly available. Hopefully, the 
work in progress on p-wave scattering will be completed, 
and our diagonalization data can shed light to the inter- 
esting problem of two dimensional dilute fermions. 

Our model of spinless fermions and hardcore bosons 
with infinite nearest-neighbor repulsion involves a signif- 
icant reduction of the size of the Hilbert space as com- 
pared to the Hubbard model. This enables us to ob- 
tain exact diagonalization results for much larger lattices 
than that can be done with the Hubbard model, and 
this also enables us to check the various analytical re- 
sults (Green function, t-matrix, diagrammatics) in the 
dilute limit with diagonalization for much larger systems 
than that has been doria.in previous works. This paper 
and a companion papcrEj on the dense limit are the first 
systematic study of the spinless fermion model in two di- 
mensions. We hope that the comprehensiveness of this 
paper can not only draw more attention to this so far 
basically overlooked model but also serve as a guide for 
diagonalization and analytical studies in the dilute limit. 
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APPENDIX A: EXACT DIAGONALIZATION 
PROGRAM 

This section describes briefly our exact diagonaliza- 
tion program. It is indebted to Refs. |6| and |7[ which 
are guides for coding exact diagonalization in one dimen- 
sion. Here we only describe the necessary considerations 
in more then one dimension, and focuses on the use trans- 
lation symmetry to reduce the problem. 

Our underlying lattice is the square lattice, and we take 
the lattice constant to be unity. The periodic boundary 
conditions are specified by two lattice vectors Ri and 
R2, such that for any lattice vector r we have r + ?iiRi + 
7121^2 = r, where ni and 712 are two integers. In Fig. 
we show two systems. The first one has Ri = (4, 0) and 
R2 = (0, 5) so the number of lattice sites is TV = 20. 
The second one has Ri = (4, 1) and R2 = (1, 5) so = 
|Ri X R2I = 19. From this example we see immediately 
the advantage of having skewed boundary conditions: we 
can have reasonably shaped systems with number of sites 
(here 19) not possible for an usual rectangular system. 
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FIG. 17: Square lattices with periodic boundary conditions: 
(4,0) X (0,5) on the left and (4, 1) x (1,5) on the right. Site 
numbers are shown, following the numbering convention, up- 
ward and rightward. 

A site order is needed to keep track of the order of the 
fermion sign. The convention that we use is starting the 
zeroth site from the lower left corner and move progres- 
sive upward and rightward following the square lattice 
structure until we encounter boundaries of the lattice de- 
fined by the periodic boundary condition vectors Ri and 
R2 (see Fig. |l^). A basis state with M particles is then 
represented by an array of the M occupied site numbers, 
with nearest neighbor s ex cluded (because V = +00 in 
our Hamiltonian Eq. (IT)). Denote such a basis state 
\n) and we have 



H\n) = -t ^ Sm\m} 



(Al) 



where Ai denotes the set of states created by hopping one 
particle in \n) to an allowed nearest-neighbor site and for 



bosons Sm = 1 always and for fermions Sm = ±1- And 
the matrix element is {m\H\n) = —Smt if m e and 0, 
otherwise. 

In order to calculate for large systems, it is necessary 
to use symmetry to block diagonalize the Hamiltonian 
matrix. In our code, we use lattice translation symmetry 
because it works for arbitrary periodic bmindaries. The 
eigenstate that we use is the Bloch stata£j 



1 



N-l 



1=0 



(A2) 



In this expression k is a wavevector (one of N , where N 
is the number of sites), R; is a lattice vector T/ is a short 
hand notation for translation by Rj, and Nn\^ is a nor- 
malization factor. The original basis states are divided 
by translation into classes and any two states in the same 
class give the same Bloch state with an overall phase fac- 
tor. What we need to do is to choose a representative 
from each class, and use this state consistently to build 
Bloch states. For a state \n) we denote its representative 
\n)- 

To compute the Hamiltonian matrix 
the Bloch states Eq. ( |A2[ ), let 

have, 



mg tne Bloch states Eq. 
representative state \n). 



us 
as 



We 

then H\Tin 



elements us- 
start fro m a 
in Eq. (^), 

- TiH\n) - 



H\n) 

~ ^mTi\m) , where we have used the fact that T/ com- 
mutes with the Hamiltonian. We have 



1=0 



(A3) 



Next because we are interested in matrix elements be- 
tween representative states, we want to connect |mk) in 
the preceding equation to |mk). If Tj(m) \ifi) = (Tj^rn) 



then Imk) = a 



(") |mk). So we have 



H\nk) = --r^ ^ Nrn-kSmCTj^m) 



(A4) 



meM 



We should note that for all to G there can be more 
than one element having the same representative |to). 
That is to say in the sum in Eq. (A4), there can be more 
than one term with |TOk). We write a new set M' — 
{to I TO G A4 and m has rep m}. Then we can write our 
matrix element equation as follows. 



(TOk|_ff |nk) 



Nn 



<^Km)e'''-'^^<-^s^. (A5) 



meM' 



Eq. ( A5) is the centerpiece of the Bloch state calculation. 
It includes many of the complications that come with 
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the Bloch basis set. (See Ref. ^ for the corresponding 
equation in one dimension.) 

Let us use M to denote the number of Bloch basis 
states for one k. M is the dimensionahty of the matrix 
that we need to diagonahze. For JV in the order of thou- 
sands, full diagonalization (with storage of the matrix) is 
done using LAPACkJ^ and a 3156 x 31Sfi matrix (7x7 
with M = 18) takes about 27 minutes.Ea For larger Af, 
the Lanczos method is implemented following the instruc- 
tions in Ref. ^ We have two options. First, we store 
information about the matrix (i.e., for each column, a 
set of (p, j{m), aj(^m)) described above that contains in- 
formation about the nonzero entries of the Hamiltonian 
matrix in this column). The M = 9 case on 7 x 7 lattice 
with = 1, 120, 744 and tolerance 10"^^ takes about 45 
minutes (32 Lanczos iterations) and uses about 1.5 GB 
of memory. This basically reaches our memory limit. 

On the other hand, we can also do Lanczos without 
storing matrix information. The same M — 9 case on 
7x7 uses only 200 MB of memory but takes more than 
four hours (263 minutes), for a larger tolerance 10"^ 
(therefore fewer Lanczos iterations, 14). Without stor- 
ing matrix information, we can calculate for larger ma- 
trices: the M = 11 case on 7 x 7, with Af = 1,906,532 
(the largest for the 7x7 system) and tolerance 10^^, 
is done in 10 hours, using less than 400 MB of mem- 
ory. The largest matrix we computed for this work is 
TV = 2,472, 147, i.e., about 2.5 million Bloch states, for 
A/ 4 on 20 X 20. This takes 10 hours and uses about 
550 MB of memory, for a tolerance of 10~^. |— , 

In addition, we have also installed ARPACKE2I that 
uses the closely related so-called Arnoldi methods and 
can obtain excited state eigenvalues and eigenvectors as 
well. If we only need information about the ground 
state, our Lanczos program is considerably faster than 
ARPACK. 

The 7x7 lattice, with maximum J\f around 2 million, 
is basically the largest lattice for which we can calculate 
eigenenergies at all fillings. The exponential growth is 
very rapid after this. The 8x8 lattice with 8 particles 
has 9, 151, 226 Bloch states, and with one more particle, 
M ~ 9, there are 30,658,325, i.e., more than 30 million 
states. 



APPENDIX B: PHYSICAL MEANING OF T{E) 

In this section we give yet another derivation of the t- 
matrix which makes more explicit the physical meaning 
of r(£:,P;q,q') Eq. (p^). 



Before we get into a lot of algebra, let us describe the 
physical idea. In scattering theory we know that the 
Born series is a perturbation series of the scattering am- 
plitude in terms of the potential. In Fig. |l^ we show 
the first three terms graphically, where the first term, 
the first Born approximation, is particularly simple-it is 
the Fourier transform of the potential. We also know 
that when the potential is weak the first few terms are 



an good approximation to the scattering amplitude, but 
when the potential is strong, we need all terms. In this 
section, we will show that our t-matrix r(i?, P,q, q') is 
the sum of all such two-body scattering terms. 



FIG. 18: The three figures represent perturbative terms in- 
volving V(q-q'), r2(£;,P;q,q') and r3(£;, P; q, q'). The 
t-matrix, r(iJ, P; q, q'), is the sum of all these terms, i.e., it 
is the sum of the ladder diagrams to infinite order. 



We start with Eq. (2.15) which we copy here for con- 
venience, 

(i?-£:(q)-£:(P-q))5(q) ^ ^J2^ici-<l')9{ci'). (Bl) 



N 
q' 

For q € Qo we break up the sum over q' into two terms 
and get. 



iE-Eo)gi<D = ^ E nq-q').9(q') 

q'eQo 

+ 4 E ^(q-q').9(q')- 



(B2) 



q'6C 



For q e Q we can rewrite Eq. (Bl) to get 



9[ 



7T5(q")- (B3) 



Plug Eq. (B3) into the second sum in Eq. (B2) and rear- 
range terms, we get. 



{E-Eo)g{ci) = ^ E nq-q')ff(q') 

q'eQo 

+ lE^2(i?,P;q,q")5(q"), (B4) 



q" 



where we have defined. 



To 



N ^ E-£{q')-£iP-<i') 
q'eQ 



k (B5) 



Now break the sum over q" in Eq. (B4) into two parts, 
and we get 

iE~Eo)giq) = ^ E nq-q')ff(q') 



q'6Qo 



+ ^ E r2(i?,P;q,q')5(q') 

q'eQo 

+ ^ E r2(£;,P;q,q').9(q')- (B6) 



q'ec 
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Plug Eq. ( p33| ) into the last term of Eq. (B6) and we get 

{E-Eo)9{<D = ^ E nq-q')ff(q') 

q'eQo 

+ ^ E r2(i?,P;q,q').9(q') 



q'eQo 



+ ^E^3(^'P'q'q')5(q'), (B7) 
q' 



where we have defined 



r3(i?,P;q,q') 



1 

7p 



E 



y(q-q")F(q"-q"')^(q" 



q",q"'eC 



{E - £(q") - £(P - q"))iE - £(q"') - £{P - q'")) 



(B8) 



Continue this process, we obtain 

{E ~ Eo)9{ci) = ^ (^q-q') 

q'eQo 

+T2{E,P;q,(^)+Ts{E,P;q,(^) + 



).9(q')-(B9) 



What we have done here is the traditional perturba- 
tion theory using iteration. Eq. ( p9| ) is the Born series 
for scattering amplitude. The first term V{q — q'), the 
Fourier transform of the potential V{r), is the first Born 
approximation. The content of higher order ter ms T 2 , T3, 
... can be obtained from their definition. Eq. (B5) say s 
that T2 involves two scatterings under V, and Eq. ( p38| ) 
says that involves three scatterings. Thus the Born 
series Eq. (|B9|).-C|an be graphically depicted at the lad- 
ders in Fig. |18|,LiJ and it involves multiple scatterings to 
all orders. Note that each term in the Born series is in- 
finite for infinite potential V. Next we will show that 



r 



summing all the terms in the series gives the t-matrix 
and the potential V cancels out, giving a finite value. 
It is easy to check that 



where Qij {E, P) is our good old Green function Eq. (p^), 

r3(£;,P;q,q') = ^ e^i-^'e-'i''^^ (^(i;, P))^, 

(Bll) 



and 



r„(i;,P;q,q') = ^ e^'i-^'e-'i'-^^ (^(-B, P))r,"'. 

(B12) 



Now plug these results into Eq. (B9), we get 



J 



{E-Eo)g{<i)= E 
q'eQo 



>^g,:q R.g-,:q' R,y 



N ^ 



(B13) 



Now we come to a formal step, 

v{i + vg + v^iGf + ...) = v{i ~ vgiE))-\ (bm) 

and the interesting result is that the infinite potential V 
cancels out, giving a finite value ^Q{E)^^. 

If we can do this formal sum, then we get from 
Eq pH), 



which is exactly our momentum space T-niatrix equa- 
tion Eg. (3.12 ) and T{E, P; q, q') is exactly our t-matrix 
Eq. (P^ . 



iE-Eo)g{ci)^ J2 r(ii;,P;q,q')ff(q') 
q'eQo 
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